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ABSTRACT 


Cloude-Pottier incoherent target decomposition (ICTD) 
and Touzi ICTD has been widely applied as a popular 
approach to interpret the scattering characteristics of a 
target in polarimetric synthetic aperture radar 
(PoISAR) data processing. However, they have a 
common drawback, i.e. proliferation of parameters 
(PoP) is unavoidable. Paladini et al. solved this 
problem by developing an orientation-invariant ICTD 
based on the coherency matrix under circular 
polarization basis. As an alternative to Paladini 
decomposition, we proposed a novel ICTD based on 
the frequently used coherency matrix under linear 
polarization basis. The proposed method can also 
avoid the problem of PoP, and avoid the ambiguity of 
alpha angle of Paladini decomposition as well. Real 
PolISAR data is processed to validate the proposed 
decomposition. 


Index Terms— Coherency matrix, ICTD, PolSAR, 
target decomposition 


1. INTRODUCTION 


Target decomposition is a popular way for target 
understanding in polarimetric synthetic aperture radar 
(PoISAR) data processing. Target decomposition 
includes two categories: coherent target 
decomposition (CTD) based on the scattering matrix 
and incoherent target decomposition (ICTD) based on 
the second order statistical coherency matrix or 
covariance matrix. As the incoherent scattering is 
general in real situation, ICTD has attracted much 
attention. ICTD originates from Huynen 
decomposition [1]. Due to its preference for symmetry 
and regularity, Huynen decomposition has not been 
widely applied [2]. Nowadays, the hotspots of ICTD 


are eigen-decomposition [3-5] and model-based 
decomposition [6-13]. 

Eigen-decomposition was pioneered by Cloude and 
Pottier [3], which was criticized by Huynen that 
serious proliferation of parameters (PoP) was 
generated [14]. Touzi also designed a successful ICTD 
based on his target scattering vector model (TSVM) 
[4], but PoP was still a problem. Paladini et al. 
proposed a novel ICTD of the coherency matrix under 
circular polarization basis, and successfully solved the 
problem of PoP [5]. A potential limitation of Paladini 
decomposition is that its alpha angle is identical to 
Cloude-Pottier alpha angle, while the helix angle is 
mixed in the alpha angle. 

For model-based decomposition, the three-component 
decomposition or four-component decomposition is 
clear in physical meaning, and easy to be implemented 
[6-9]. However, their disadvantage is information loss. 
Chen et al. proposed a novel model-based 
decomposition with separate orientation angle for odd- 
and double-bounce models [10], but we think the helix 
component should be modeled as the asymmetry of 
odd- and even-bounce scatterers, not as an 
independent scattering component. Van Zyl et al., Cui 
et al., and Wang et al. proposed nonnegative eigen- 
decomposition [11-13], but the interpretation of 
scattering mechanisms is partially based on eigen- 
decomposition, not on scattering models. 

As an alternative to Paladini decomposition, our ICTD 
which bases on the coherency matrix under linear 
polarization basis also avoids PoP, and also avoids the 
alpha angle ambiguity of Paladini ICTD. 


2. MATRIX DEFINITIONS 


For coherent scattering, the full polarimetric scattering 
information is contained in the scattering matrix 


s-| $ a (1) 


where the subscript HV denotes vertical polarization 
transmission and horizontal polarization reception. 
Two of frequently used scattering vectors are the 
scattering vector kp in Pauli basis and the scattering 
vector kc in circular polarization basis 


1 
kp = [Sim + Sw Siu = Sy Bwl CO) 


t 
ke =| Six V251, Sr, | (3) 
where the superscript t denotes matrix transposition, 
and the subscripts L and Lı denote left circular 
polarization and left orthogonal polarization. 
For incoherent scattering, the full polarimetric 
information is contained in the second order statistical 
coherency matrix T in linear polarization basis or 
coherency matrix G in circular polarization basis 


T= (kp Kh) (4) 


G=(ke-ke) (5) 
where <-> denotes time average or spatial average, 
and the superscript + denotes complex conjugation and 
transposition. Both T and G have nine degrees of 
freedom (DoFs). 

The Cloude-Pottier, Touzi, and the proposed ICTDs 
are based on the coherency matrix T, while the 
Paladini ICTD is based on the coherency matrix G 


3. THE PROPOSED ICTD 


The proposed ICTD is based on eigen-decomposition 
of the coherency matrix T 


T=UAU! (6) 
where A is the diagonal matrix whose diagonal 
elements are the eigenvalues 2,, 42, and 43, with 
A2423 assumed, and U is comprised of three 
columns of eigenvectors kı, k2, and k3. 

The matrix U is modelled by the multiplication of six 
unitary transformation matrices, and the DoF of each 
unitary transformation matrix is one 
6 
U=[k ky kgl=[[SUGy,) (7) 
i=l 
In the following subsections we will detail the models 
of the matrices A and U. 


3.1. Information extraction from the matrix A 


The matrix A has the form A=diag(A), 42, 43). It has 
three independent parameters. For convenient physical 
information extraction, other three physical 
meaningful parameters are here suggested to replace 
the three eigenvalues, i.e. the total power SPAN, the 
scattering entropy H, and the anisotropy A 


SPAN = 4 +% +4, (8) 
3 A 
H =—)_p;log3 p; , with p; =———*——__ (9) 
a i ' At’ +, 
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The SPAN is the total power backscattered from the 
target, which was introduced by Cao ef al. into the 
classification applications [15]. The entropy H is an 
indicator of the scattering randomness. The anisotropy 
A is also physical meaningful. 


3.2. Model of the matrix U 


As the matrix A has three DoFs, the matrix U has only 
six DoFs. 


3.2.1. Model of the dominant eigenvector kı 

For three dimensional complex vector k,, there are six 
DoFs. Two DoFs are lost because of 1) absolute phase 
indetermination and 2) unitary restriction [5]. Thus 
four independent parameters are used to model ky. 

We model k; using Touzi’s TSVM as 


1 0 0 COS A, COS 2T,, 
ky =|0 cos2y —sin2y sin a e)?as 
0 sin2yw cos2y || —jcosa, sin2r,, 
1 0 0 1 0 0 
=|0 cos2y —sin2y || 0 elas 0) 
O sin2y cos2y || 0 0 1 
cos2t, 0 —jsin2z,, ||cosa, —sina, O]| 1 
0 1 0 sina, cosa, 0O/||0 
—jsin2t,, 0  cos2z,, 0 0 110 
1 
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0 
(10) 


The parameter y is the orientation angle of k,. The 
parameter 7,,is the helix angle of kı, which measures 


the asymmetry of the target represented by kı. The 
parameter a, is the scattering type magnitude, which is 
an important parameter used to discriminate the odd- 
and the even-bounce scatterings. The parameter ®,, is 
the scattering type phase, which can be used to 
discriminate dipoles and quarter waves and has been 
used for ship recognition and other applications [4]. 
The alpha angle of Paladini ICTD is identical to the 
Cloude-Pottier alpha angle. But the Cloude-Pottier 
alpha angle is a mixer of alpha angle and helix angle, 
it is ambiguous in physical meaning [4]. The 
parameter a, of TSVM or our ICTD can avoid such 
ambiguity. 

The four parameters have clear physical meanings. 
They eliminate some ambiguities of the Cloude-Pottier 
alpha-beta model [4]. 


3.2.1. Model of the subdominant eigenvectors 
As the matrix U has six DoFs, and the dominant 
eigenvector has four DoFs. Thus only two DoFs are 
left in the two subdominant eigenvectors k and kz. 
After conducting the following unitary 
transformations, we can get [5] 

X = SU (~a, )SU (~T m )SU (~® gs )SU (-w)U 


1 0 0 
(11) 
= (0) x X2 P 
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The subdominant eigenvectors are modelled by two 
unitary transformation matrices as follows 


1 0 0 
0 x x» |=SUDSU(Y) 
O x3 X4 
1 0 0 1 0 0 (12) 
=|0 eT o |10 cos2¥ -sin2¥ |. 
0 0 elI |0 sin2¥ cos2¥ 


The parameter ¥ measures the relative magnitude of 
the second and third elements of the eigenvector, 
which is interpreted as the local orientation angle of 
the target. The parameter I introduces some phase 
difference containing the information of the local helix 
of the target. The physical basis of these two 
parameters needs further validation. 

In [16], the last two parameters are interpreted as the 
“relative orientation” and “relative helicity”, and they 
are both depolarization parameters. 


In low entropy scattering situation, the dominant 
eigenvector described by four parameters accounts for 
the dominant or global scattering characteristics. The 
subdominant eigenvectors described by the remained 
two parameters account for the local scattering 
characteristics of the target. In high entropy scattering 
situation, the scattering power of the subdominant 
eigenvectors are comparable to that of the dominant 
eigenvector, and more information is contained in the 
last two parameters compared to the low entropy 
situation. 

By the proposed ICTD, nine independent parameters 
are decomposed from the coherency matrix, i.e. SPAN, 
H, A, Y, Tm, Os, Pas, Y, and I. Thus the proposed ICTD 
can avoid the problem of PoP. All the decomposed 
parameters have explicit physical meanings. 


4. EXPERIMENTAL RESULTS 


In this section, RADARSAT-2 C-band PolSAR data of 
San Francisco Bay area are processed to demonstrate 
the inherent physical meanings of the nine 
decomposed parameters. 

The data are in single-look complex format. A refined 
Lee filter is then applied to suppress the speckle and 
get the coherency matrix format data. Only a part of 
the scene is utilized. 

The total power SPAN is shown in Fig. 1(a), as one 
can see, dihedrals with specular scattering usually 
have higher backscattered powers. The scattering 
entropy H is shown in Fig. 1(b), from which one can 
see that the park area and the oriented urban area with 
higher scattering randomness have higher entropy. The 
anisotropy A is shown in Fig. 1(c). For ocean area and 
urban area, the minimum eigenvalue is very small, and 
thus A is large. The orientation w of the dominant 
eigenvector is shown in Fig. 1(d), where the blue color 
area reasonably corresponds to the oriented urban area. 
The helix angle tm of the dominant eigenvector is 
shown in Fig. 1(e), where small value of helix angle is 
exhibited for the whole scene. The scattering type 
magnitude a, is shown in Fig. 1(f), which is close to 
m/2 for urban area, while close to zero for ocean area. 
The scattering type phase ®,, is shown in Fig. 1(g), 
which also contains physical information as addressed 
before. The local orientation ¥ is shown in Fig. 1(h), 
which is large for oriented urban area and park area 
because the sub-scatterers in the cell have diverse 
orientations. The local helix angle I is shown in Fig. 
1i), and it is noise generally. 
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5. DISCUSSION 


Similarities and differences of the Paladini ICTD and 
the proposed ICTD are addressed. 

First, the similarities are 1) they both avoid the 
problem of PoP; 2) they are both serial (multiplicative) 
decompositions, while other common decompositions 
which are parallel (additive) decompositions. 3) the 
helix angles of both ICTDs have ambiguities. The 
helix angle cannot describe the asymmetry when a, 
equals 2/2 for proposed ICTD, while the helix angle 
cannot describe the asymmetry when alpha angle 
equals zero for Paladini ICTD. 

Second, the differences are 1) our proposed ICTD 
bases on the frequently used coherency matrix T, 
while Paladini ICTD bases on the coherency matrix G; 
2) the dominant eigenvector model for our ICTD is 
based on TSVM, and it eliminates the alpha angle 
ambiguity of Paladini ICTD; 


6. CONCLUSION 


A novel serial (multiplicative) ICTD of the coherency 
matrix is proposed with nine decomposed independent 
parameters obtained and all of them have explicit 
physical meanings. The proposed decomposition 
successfully avoids the problem of PoP. Real PolSAR 
data are processed to demonstrate the effectiveness of 
the proposed approach. The classification method 
based on these nine parameters is under study. 
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Figure. 1. Results of the parameters (a) SPAN, (b) H, (c) 
A, (d) Y, (e) Tm, (f) Os, (g) Das, (h) Y; and (i) T. 


